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ABSTRACT 

We present calculations indicating the possibility of a new class of Type I supernovae. 
In this new paradigm relativistic terms enhance the self gravity of a carbon-oxygen 
white dwarf as it passes or orbits near a black hole. This relativistic compression can 
cause the central density to exceed the threshold for pycnonuclear reactions so that 
a thermonuclear runaway ensues. We consider three possible environments: 1) white 
dwarfs orbiting a low-mass (~ 10 — 20 M ) black hole; 2) white dwarfs encountering a 
massive (~ 1 — 3 x 10 3 M ) black hole in a dense globular cluster; and 3) white dwarfs 
passing a supermassive (~ 10 6 — 10 9 M ) black hole in a dense galactic core. We estimate 
the rate at which such events could occur out to a redshift of z ~ 1. Event rates are 
estimated to be significantly less than the rate of normal Type la supernovae for all three 
classes. Nevertheless, such events may be frequent enough to warrant a search for this 
new class of supernova. We propose several observable signatures which might be used 
to identify this type of event and speculate that such an event might have produced the 
observed "mixed-morphology" Sgr A East supernova remnant in the Galactic core. 
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1. Introduction 

Type I supcrnovae are currently objects of intense 
interest. The observed luminosity of Type la super- 
novae (SNIa) out to high redshift offers some of the 
most convincing evidence for the acceleration of the 
universe due to an unknown dark energy (Garnavich 
ct al 1998; Riess et al. 1998; Perlmutter et al. 1998). 
The distance redshift relation for the SNIa standard 
candle is also an important component in the deter- 
mination of the present cosmic expansion rate (Freed- 
man et al. 2001). As such, it is important to carefully 
scrutinize any possible variations in the paradigm for 
SNIa explosions which might occur at large redshift. 
In this paper we explore one such possible new mech- 
anism for SNIa's. 

Type I supernovae are characterized by the absence 
of hydrogen lines in their spectrum. Type la super- 
novae, in particular, are believed to result from the 
accretion of material onto the surface of a carbon- 
oxygen white dwarf. If the accretion rate is suffi- 
ciently rapid, mass accumulates until the star ap- 
proaches the Chandrasekhar limit. As the interior 
density exceeds the threshold for pycnonuclear reac- 
tions, an off-center runaway deflagration/detonation 
wave of thermal carbon burning ignites under strongly 
degenerate conditions. Since the degeneracy pressure 
support of the core is almost independent of temper- 
ature, the heating of the core does not cause the star 
to expand and cool. Therefore, carbon burning rates 
increase until a thermonuclear explosion eventually 
engulfs a significant fraction of the star. The key to 
this paradigm is that the density at some time ex- 
ceeds the threshold for thermonuclear ignition. For 
typical Type la models the ignition density ranges 
from pi gn 2 — 4 x 10 9 g cm -3 (at a temperature 
T~ 5 x 10 8 K) with the lower range of ignition den- 
sity favored by nucleosynthesis arguments (Iwamoto 
ct al. 1999). 

In this paper we propose a new mechanism by 
which the ignition density of a white dwarf can be 
achieved. We show that relativistic enhancements 
of the self gravity of a white dwarf occur whenever 
the star is accelerating in a background gravitational 
field. As the velocity with respect to the background 
field approaches a significant fraction of the speed of 
light (as will occur near a black hole) these relativis- 
tic effects can cause the central density to exceed the 
threshold for nuclear ignition even for white dwarfs 
which arc well below the Chandrasekhar mass. Alter- 



natively, the ignition densities might be achieved by 
the excitation of tidally-induced resonance/disruption 
(Khokhlov & Melia 1996; Rathore & Blanford 2000) 
within the white dwarf, or tidal compression by the 
Carter-Luminet effect (Carter & Luminet 1982; La- 
guna, Miller & Aurek 1994). In any of these sce- 
narios, the result is the same, i.e. the ignition of a 
thermonuclear explosion in a deep gravitational well. 
The focus of the present work, however, is on the rel- 
ativistic compression effect (Wilson & Mathews 2000; 
Wilson 2002). 

We note that other authors have modeled hydrody- 
namic effects which occur as stars encounter a black 
hole. For example, Fryer et al. (1999) have consid- 
ered the tidal disruption of a white dwarf near a low 
mass black hole as a source of gamma-ray bursts. 
Their simulations were based upon Newtonian hy- 
drodynamics which would not incorporate the gen- 
eral relativistic effect described here. Similarly, Ayal, 
Livio, and Piran (2000) have simulated the tidal dis- 
ruption of a normal solar type main-sequence star 
near a massive (m <~ 1000 M©) black hole. Extended 
main sequence stars indeed should tidally disrupt be- 
fore the relativistic effects described here can induce 
compression. The effect described here is unique to 
compact objects such as white dwarfs or neutron stars 
and is only apparent in a fully relativistic simulation. 

We consider three possible environments in which 
such events might occur. These are: 1) white dwarfs 
orbiting a low-mass black hole (hereafter LMBHSN); 
2) white dwarfs encountering a massive (M > 1000 
M© black hole in a globular cluster (hereafter MBHSN) ; 
and 3) white dwarfs passing a supermassive (M > 
10 6 M© black hole in a dense galactic core (here- 
after SMBHSN). In what follows we utilize general 
relativistic hydrodynamic numerical simulations de- 
scribed in Sections 2 and 3 to determine the condi- 
tions under which the central density can be made to 
exceed the ignition threshold for these three paradigms. 
In Section 4 we then make a rough analysis of the 
event rates at which these various scenarios might oc- 
cur. We conclude with some discussion as to how such 
events might be identified observationally. 

2. White-Dwarf Compression 

It has been noted in studies of colliding (Wilson 
2002) and orbiting (Mathews & Wilson 2000) neu- 
tron stars that general relativistic effects can over- 
whelm the stabilizing tidal forces and cause the cen- 



2 



tral density of compact objects to increase. This phe- 
nomenon was shown (Mathews, Marronetti & Wilson 
1998; Mathews & Wilson 2000) to scale with the mag- 
nitude of the spatial component of the four velocity 
of the star relative to the background field for condi- 
tions in which there exists no Killing vector by which 
the three-velocity can be made to vanish. 

Here we point out that these conditions are also 
met for white dwarfs orbiting or passing near a black 
hole. Thus, a similar scaling of interior density with 
four velocity should be expected for the case of white 
dwarfs accelerating under the influence of the grav- 
itation potential of a massive black hole. We show 
by numerical calculations that the increase of central 
density is sufficient to induce a thermonuclear run- 
away and lead to a new class of Type I supernovae. 
These supernovae are unique in that they would most 
likely be found in association with a black hole. The 
case of MBHSNc or SMBHSNc would also be asso- 
ciated with dense globular clusters or galactic cores, 
respectively. They would also be distinguished from 
the usual Type la event in that they may involve a less 
massive white dwarf and a rapid ignition/detonation 
in their centers. 

2.1. The Model 

We analyze the induced compression in the con- 
text of a general relativistic treatment of the hydro- 
dynamic and field evolution of a white dwarf in the 
vicinity of a black hole. Although an orbiting system 
is inherently nonaxisymmetric, the essential features 
of this system can be modeled (with slight modifica- 
tions) while imposing axisymmetry. The implemen- 
tation of axisymmetry does not lead to significant in- 
accuracies, but it allows us to rapidly explore a broad 
range of parameter space for white-dwarf plus black- 
hole systems as follows. 

An exact axisymmetric general relativistic hydro- 
dynamic model has been previously developed by Wil- 
son (1979; 2002) (see also Hawley, Smarr & Wilson 
1984) to describe the head-on collision of two neu- 
tron stars. Here we have modified this model to treat 
the essential features of a white-dwarf plus black-hole 
system. 

We start with the ADM formalism (Arnowitt, 
Deser & Misner 1962) whereby, the metric is split 
into space-like and time-like components 

ds 2 = -{a 2 - (3 t [3 l )dt 2 + 2(3 l dx l dt + 7 y -cfeW , (1) 

where a is the lapse function, p l is the shift vector, 



and 7y is the metric for the three-space (we use ge- 
ometrized units. G = c = 1). For the present pur- 
poses, we are concerned with either quasi-stationary 
orbits or passing orbits near the point of closest ap- 
proach. In either case, the P r and f3 z metric compo- 
nents are very small for the systems of interest here 
and can be neglected in what follows. In a full nonax- 
isymmetric treatment, the dominant contribution to 
the shift vector will be from nonaxisymmetric motion 
about the rotation axis of the orbit plane. Neverthe- 
less, an axisymmetric model can still be maintained 
by modified hydrodynamic equations as described be- 
low. 

In essence there are two effects of orbital motion 
on the shift vector: In a comoving gauge (Wilson & 
Mathews 2003) the dominant contribution is from an 
overall rotation of the coordinates. The explicit inclu- 
sion of this term naturally generates the appropriate 
centrifugal and Corriolis forces associated with a ro- 
tating coordinate system. There is also a small Lense- 
Thirring frame-drag correction (Wilson, Mathews & 
Marronetti 1996) for the relativistic orbit precession. 
The frame-drag component can be deduced from an 
equation of the form, \7 2 f3 drag sa Wwa^WpwDU 
(Wilson & Mathews 2003). For the systems of inter- 
est here, the source term is small so that frame-drag 
is only a small correction to the orbital motion. This 
is consistent with numerical results for neutron-star 
binaries near their last stable orbit (Mathews & Wil- 
son 2000). Hence, for our purposes the frame-drag 
can be neglected. 

The coordinate rotation term can be eliminated 
by choosing an instantaneously nonrotating coordi- 
nate system at the orbit periapse and adding effective 
centripetal force terms to the equations of motion as 
described below. 

Thus, for this application we adopt a quasi-Eulerian 
(J3i = 0) gauge (Wilson & Mathews 2003). We further 
adopt a conformally flat metric whereby 7^ = <fi 4 5ij. 
This has been shown (Cook et al. 1993; Mathews & 
Wilson 2000; Wilson 2002) to be sufficiently accurate 
for our purpose. 

Thus, in cylindrical coordinates we write an effec- 
tive metric , 

ds 2 = -adt 2 + 4> 4 [dr 2 + dz 2 + r 2 d6 2 } , (2) 

where in what follows we take the z coordinate to 
be the distance along the line of centers between the 
black hole and the white dwarf. The r coordinate is 
the radius perpendicular to the z axis. This metric 
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leads (Wilson, Mathews & Marronetti 1996; Wilson 
2002; Wilson & Mathews 2003) to simple Poisson- 
like elliptic constraint equations which can be solved 
to find the the conformal factor and lapse function. 
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where the source term includes Lorentz contracted 
state variables, D = pW, E — epW, with p the local 
proper mass density, e the internal energy per unit 
mass, and W a Lorentz-like factor, W = \/l + U l Ui 7 
where U l denotes spatial components of the four ve- 
locity. 

The quantity T = 1 + P/ep is an equation of state 
index. We use a standard white dwarf equation of 
state (e.g. Salpcter 1961; Hamada & Johnson 1961) 
to specify e and/or T as a function of density. For sim- 
plicity we only retain the dominant degenerate elec- 
tron part. 

In equations (3) and (4) the extrinsic curvature 
KijK" terms result (Wilson & Mathews 2003) from 
derivatives of the relativistic frame drag contained 
in the /?'. In close neutron-star binary systems it 
has been shown (Mathews & Wilson 2000) that this 
term contributes insignificantly to the metric source 
terms, (/ dVK i3 K^ /8tt ~ 10- 5 M Grav ). In what fol- 
lows these extrinsic curvature terms can therefore be 
neglected. Including these terms would only slightly 
enhance the compression effect described below. 

2.2. Axisymmetric Hydrodynamics 

To model a white dwarf near a black hole, a cen- 
trifugal term is added to the axisymmetric momen- 
tum equations to keep the stars apart. Thus, the 
hydrodynamic equations are slightly modified from 
those presented in (Wilson 2002; Wilson & Mathews 
2003) The momentum equation becomes, 
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where a = D + YE, and the covariant coordinate 
momentum density is Si = <jJJi, and U 2 = UiU 1 . 

The last term in Eq. (5) is a centrifugal accelera- 
tion. It accounts for the conserved motion along an 
unspecified third coordinate perpendicular to z. Ob- 
viously, setting this term to zero recovers the simple 
head on collision. If this term is sufficiently large, it 
can stabilize the stars in orbit, or it can limit motion 
along z to a distance of closest approach. 

It is to be noted that this component must be in- 
cluded in the UiU 1 summations above and in the met- 
ric fields. For example, at closest approach, U 2 = 
U±U ± . Thus, the inclusion of this term increases the 
source terms [cf. Eqs. (3) and (4)] for the solution of 
the metric fields. 

The remaining hydrodynamic equations assume 
their usual form. The continuity equation is, 
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The equation of internal energy evolution is, 
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2.3. Enhanced Self Gravity 



The term with aad In a/dx l in Eq. (5) is the 
the relativistic analog of the Newtonian gravitational 
force. The general relativistic condition of hydrostatic 
equilibrium for each star can be inferred by taking the 
time-stationary limit of the momentum equation (5), 
and transforming to instantaneously comoving three 
space coordinates with a grid velocity x l g 
such that the three velocity V j 
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0. This gives 
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Note, that in the weak-field limit, <\> 2 
reduces Eq. (8) to 



a which 



dSj, ao dlncj) 1 8 

dt + bbl dt + dfidxi 



>SiV j + a 



dP 

dx % 



dP 

dx l 



d In a 

r — 

dx l 



1 + 2C/' 



(9) 



4 



from which it is manifestly apparent that the terms 
in the square brackets of Eq. (8) do not cancel. The 
enhanced self gravity derives from the additional U 2 
dependent terms of Eqs. (8) and (9). 

3. Model Calculations 

As noted above, we consider several possibilities in 
this study. One is a white dwarf in a bound orbit 
around a low-mass black hole. This situation could 
result, for example, from the evolution of a binary in 
which one member was a massive progenitor which 
formed a black hole. Alternatively, a white dwarf 
might be captured by a black hole through three-body 
interactions in a dense cluster. In either case, the 
white dwarf orbit would decay by gravitational radia- 
tion until it approaches high orbital velocity. Another 
possibility is that of a white dwarf randomly passing 
a massive >1000 M black hole in a dense globular 
cluster or galactic core. 

To start each calculation an initial distribution of 
density and internal energy is taken (e.g. a white 
dwarf in hydrostatic equilibrium) and with initial ve- 
locities corresponding to either nonradial infall from 
infinity, or of an amount required to stabilize the or- 
bit, i.e. an angular momentum is chosen to present a 
centrifugal barrier. In the case of a simple bound cir- 
cular orbit, this term was simply adjusted to keep the 
white dwarf at a fixed distance from the black hole. 
For the case of passing orbits, this term was chosen 
to produce a particular distance of closest approach. 
The stars were then evolved and the central density 
monitored until the white dwarf was either tidally dis- 
rupted, attained the critical density for thermonuclear 
burn, or passed the distance of closest approach. 

For our purposes, a thermonuclear runaway is con- 
servatively expected to initiate once the density ex- 
ceeds 3 x 10 9 g cm' 3 . Table 1 summarizes key prop- 
erties of white dwarfs of various mass in bound cir- 
cular orbits around low-mass black holes. Tables 2 
and 3 summarizes some calculations of passing orbits 
of massive and supermassive black holes, respectively, 
at the points of closest approach. For illustration, col- 
umn 3 in these tables gives the central density for 
isolated white dwarfs of the indicated mass. 

Most of these table entries have been selected be- 
cause they represent systems which have just achieved 
the ignition density in their centers at the distance 
of closest approach. The two entries of a 1.2 M Q 
white dwarf orbiting a low mass black hole, however, 



are just below our adopted ignition density. All of 
the stars in Table 1 have large tidal distortion even 
though their central density is quite high. 

Figures 1 to 3 illustrate the competition between 
tidal distortion and relativistic compression. These 
figures arc for a typical 0.6 M Q white dwarf passing 
near various massive black holes. Note, that the ra- 
dius of an isolated 0.6 M Q white dwarf is ~ 9000 km. 
Hence, these stars are significantly compressed. The 
black hole resides at <~ 10 4 km below these figures 
as noted in the captions. All three cases have been 
chosen because the white dwarf just achieves the ig- 
nition threshold in its center near the point of closest 
approach to the black hole. These therefore represent 
the maximum impact parameter for which explosion 
can occur. 

Figure 1 is for a white dwarf near a 1000 M© black 
hole. In this case, the white dwarf is substantially 
tidally distorted. Nevertheless it has compressed (by 
more than a factor of 1000 in density and a factor of 
10 in radius) so that the center reaches the ignition 
threshold. Figure 2 is for a 1500 M Q black hole. It 
shows less tidal distortion since it reaches the central 
ignition density farther from the star. Similarly, Fig- 
ure 3 for a white dwarf even farther from a 3000 M Q 
black hole exhibits almost no tidal distortion at all. 
These figures clearly illustrate that for such systems, 
the tidal forces (which tend to decrease the central 
density) are not strong enough to prevent an increase 
in central density as the white dwarf approaches the 
black hole. 

From these simulations the following conclusions 
can be made regarding this paradigm: 

• Substantial increases in central density can oc- 
cur as a star passes close to a black hole. For 
example, the central density of a 0.6 M Q white 
dwarf in isolation is 3.2 x 10 6 g cm~ 3 . This den- 
sity is increased by nearly a factor of 10 3 when 
passing a massive black hole. 

• The ignition density can be achieved for smaller 
velocities for white dwarfs that are closer to the 
Chandrasckhar mass. 

• For low-mass black holes, pycnonuclear igni- 
tion requires the orbit to be relatively close 
(~ few x 10 3 km), i.e comparable to the radius 
of the white dwarf. Also, the white dwarf must 
be in the high end of the observed mass distri- 
bution, Mwd > 1-2 M Q . However, For a black 
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Table 1 

White Dwarfs Orbiting a Black Hole 



M WD (M ) M BH (M ) pf (10 9 g cm- 3 ) p c (10 9 g cm" 3 ) z min (km) 



U 2 



1.20 
1.20 
1.30 
1.30 



10 
20 
10 
20 



0.12 
0.12 
0.37 
0.37 



1.2 
1.6 
2.9 
3.0 



4.82 x 10 3 0.0031 

6.52 x 10 3 0.0045 

4.21 x 10 3 0.0035 

5.77 x 10 3 0.0051 



Table 2 

White Dwarfs Passing a Massive Black Hole 



M WD (M ) 


M BH (M ) 


p™ (10 9 g cm" 3 ) 


p c (10 9 g cm" 3 ) 


z mm (km) 


U 2 




0.60 


1000 


0.0032 


3.4 


0.120 x 10 b 


0.248 


0.60 


1500 


0.0032 


3.2 


0.193 x 10 5 


0.231 


0.60 


3000 


0.0032 


3.9 


0.393 x 10 5 


0.226 


0.80 


3000 


0.0091 


3.6 


0.706 x 10 5 


0.126 


1.00 


3000 


0.029 


3.1 


1.24 x 10 5 


0.0715 


1.20 


3000 


0.12 


3.1 


2.81 x 10 5 


0.0316 


1.30 


3000 


0.37 


3.2 


5.50 x 10 5 


0.0160 


1.40 


3000 


1.9 


3.2 


33.0 x 10 5 


0.0027 



Table 3 

White Dwarfs Passing a Supermassive Black Hole 



M WD (M ) M BH (M ) pf (10 9 g cm- 3 ) p c (10 9 g cm" 3 ) z mm (km) U 2 



0.60 10 9 0.0032 3.06 1.5 x 10 10 0.191 

0.70 10 9 0.0056 3.29 2.0 x 10 10 0.151 

0.80 10 9 0.0091 3.05 2.5 x 10 10 0.118 

0.90 10 9 0.010 2.95 3.3 x 10 10 0.0905 

1.00 10 9 0.029 3.13 4.3 x 10 10 0.0686 
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hole of rriBH > 1000 M Q sufficient compression 
can occur even for a 0.6 M© white dwarf, and 
ignition can occur much further from the black 
hole where tidal effects are minimal (cf. Fig. 3). 

Hence, we conclude that the most likely environ- 
ment for this effect is in a dense galactic core or glob- 
ular cluster in which there is a good chance that a 
white dwarf can pass relatively close to a massive or 
supermassive black hole. 

4. Event Rates 

4.1. Normal Type la Supernovae 

It is important to have an estimation of the rates 
of such new events relative to those of normal SNIa's. 
The observable event rate for normal SNIa explosions 
out to a rcdshift of z w 1 is rather significant 

Nsni = N gal R Ia « 10 9 x 10- 2 « 10 7 yr- 1 , (10) 

where N ga i is the number of galaxies out to z w 1 , and 
Rj a is the average rate per galaxy at which carbon- 
oxygen white dwarfs accrete enough material to ex- 
ceed the critical interior density for the ignition of an 
explosion. This rate is observed to be about 1 per 
century per galaxy (Kobayashi et al. 2000). It is in- 
structive, however, to review how one estimates this 
rate theoretically: 

If we adopt the current view that Type la super- 
novae arise from the binary accretion of material onto 
a white dwarf. There are two scenarios which have 
been debated in the literature (Nomoto et al. 1994; 
Kobayashi et al 2000). The most widely adopted from 
a theoretical modeling point of view (and the one 
that will be adopted here) is that of a single degen- 
erate white dwarf which accretes hydrogen via mass 
transfer from a normal red giant or main sequence 
companion. If the mass transfer is sufficiently rapid, 
the interior density rises to the ignition density before 
material can be ejected via normal nova eruptions. A 
second possibility (Iben & Tutukov 1984; 1985) is that 
of a doubly degenerate carbon-oxygen white-dwarf bi- 
nary. The eventual merger of the two white dwarfs 
leads to a single system of order the Chandrasekhar 
mass which can ignite and explode. 

The relative rates of cither of these paradigms can 
be deduced from the relative probability of having a 
close interacting binary. For this we utilize a slightly 
modified version (Mathews, Bazan & Cowan 1992) 



of an expression first proposed by Greggio & Renzini 
(1983). 

Ria(t) = / dM B <t>{M B ) / dqf(q) 

J Mi Jqi 
j-a h (t-T,M B ,q) 

x / daP(a)ip(t - t) , (11) 

J ai (t,MB,q) 

where 4>(Mb) is the initial mass function (IMF) for 
a binary of total mass, Mb, f{q) is the distribution 
function for the ratio of the secondary to primary 
mass, q = m 2 /mi, P(a) is the distribution of ini- 
tial separations, a, and tp(t — r) is the star formation 
rate corrected for the evolutionary time r for the sec- 
ondary star to evolve to a white dwarf and for the 
binary to evolve to a supernova. As in (Mathews, 
Bazan & Cowan 1992) we take f(q) oc q x with x < 2, 
since distributions with exponents > 2 give incorrect 
ratios of single lined to double lined spectroscopic bi- 
nary systems (Kraicheva et al. 1979). This favors 
mass-symmetric systems. The initial separation dis- 
tribution is taken by similar arguments (Tutukov & 
Yungelson 1979) to be P(a) oc a -1 , which favors close 
systems. 

The upper and lower limits to the mass integral 
are simply the largest and smallest total binary mass, 
respectively, allowed by the scenario and evolutionary 
time scale. For a single degenerate white dwarf SNIa, 
one typically has, 1 < m2 < 3 M Q and 3 < mi < 8 
M as the carbon-oxygen white dwarf progenitor. 
This implies, Mi = 4, M h = 11 and q t = 1/8, q h = 1. 
The limits on the separation distance are determined 
by the presumed scenario for Type la supernovae as 
discussed above. For our purposes, the lower limit 
a t <~ 10 AU is determined (Mathews et al. 1992) by 
the requirement that the stars not be too close during 
the common-envelope asymptotic-giant-branch phase 
(Iben & Tutukov 1984; 1985) of the white dwarf pro- 
genitor to prevent the formation of a CO white dwarf. 

4.2. White Dwarfs Orbiting a Black Hole 

One might think that a white dwarf orbiting a com- 
mon low-mass black hole could be the most likely can- 
didate for this event. After all, a number of candidate 
binary systems have been identified (Cowley 1992) 
which probably contain a black holes as one member. 
Indeed, mass transfer from a close companion star 
is the primary means to identify a black hole candi- 
date. Hence, it seems quite likely that a number of 
close binary systems exist in which the higher-mass 
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progenitor star has evolved to a black hole while the 
lower-mass member is now a white dwarf. In this case, 
the inspiraling of the white-dwarf companion is an in- 
evitable outcome and the possibility of a compression- 
induced explosion needs to be considered. Here, we 
estimate the rate of such events relative to normal 
SNIa events. We find that this is a possible, though 
rare, phenomenon. 

The relative rate of exploding white-dwarf, black- 
hole binaries compared to that of normal Type I su- 
pernovae can be deduced from Eq. (11). Ignoring 
the logarithmic dependence on relative separation, 
and assuming instantaneous recycling along with a 
Salpeter IMF, the ratio of rates becomes: 



R 



LMBHSN 
RsNla 



-1.3 



M h 1 ' ? ']lmbhsn 



[Mf 1 - 3 



x+1 
max 



, n \LMBHSN 



[Umax IminlSNla 



(12) 



Adopting the parameters above for the carbon-oxygen 
white dwarf progenitor, and progenitor masses for a 
low-mass black hole in the range of 20 < mi < 60, 
we have Mi = 23, M h = 68, q min = 0.04, and q max = 
0.40 for the LMBHSN scenario. Setting x = 1 then 
implies a most optimistic rate of 



Rlmbhs 



N 
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10" 



(13) 



Requiring that the white dwarf have a mass M\y d > 
1.2 M and that the progenitor star produce a black 
hole with M BH > 20 M Q (as indicated in Table 1) 
will decrease this rate by at least another order of 
magnitude. Hence, we conclude that this type of su- 
pernova is probably much less frequent than normal 
Type la supernovae. The main reason is the simple 
fact that the high-mass progenitors of a black hole 
are more rare than the low-mass progenitors which 
can produce SNIa's. Also, the required massive white 
dwarfs are rare, and the required binary systems with 
a high mass asymmetry are less common. 

5. White Dwarfs Passing a BH 

The other plausible candidate environments for the 
explosion scenario considered here are those of white 
dwarfs on passing orbits near a massive (itibh ~ 10 3 
Mq) black hole in the center of a globular cluster or 
a supermassive (itibh ~ 10 6 — 10 9 M Q ) black hole 
in a dense galactic core. We show here that is possi- 
ble for event rate in these environments to be greater 



than the event rate due to low-mass black-hole/ white- 
dwarf binaries. Supermassive black holes, for exam- 
ple, appear to be present in every galaxy. They also 
tend to be found in high density stellar environments 
in which a large encounter rate is possible. 

5.1. Massive Black Hole in Globular Clusters 

At the present time there is no firm observational 
evidence for massive black holes in globular clusters 
(GCs) even though much effort has been expended 
searching for them. Although there have been re- 
cent claims (e.g. Gebhardt, Rich & Ho 2002; van der 
Marel et al. 2002) for their detection, it has since been 
demonstrated (cf. Baumgardt et al. 2003; Ferraro et 
al. 2003; McNamara 2003) that conventional inter- 
pretations (without a massive black hole) can match 
the observations equally well. Hence there is no firm 
empirical basis to conclude that most globular clus- 
ters harbor an intermediate mass black hole. 

Nevertheless, there are plausible theoretical argu- 
ments as to why such objects should exist in at least a 
fraction of globular clusters (Miller & Hamilton 2002; 
Zheng 2001). For example, a strong correlation has 
been noted (Gebhardt et al. 2000; Ferrarese & Mer- 
ritt 2000) between core black-hole masses and the 
dispersion velocity in the central bulge of the galax- 
ies. Observed black-hole masses scale (Tremaine et 
al. 2002) as, 



M, 



BH 



TO 



M P 



(14) 



where a = 4.02 ± 0.32 and M w 1.35 x 10 8 for 
(T = 200 km s _1 . It has been shown (Adams et al. 
2003) that this scaling simply follows because collapse 
to the black-hole mass is halted once the centrifugal 
radius of the collapse flow exceeds the capture radius 
of the central black hole. A similar collapse should 
account for formation of any spherical stellar system. 
Hence, it is reasonable to extrapolate this relation 
to globular clusters. For the 22 Milky- Way globu- 
lar clusters considered in (Zheng 2001) ~ 10% would 
have implied black hole masses of the order of 10 3 
M , but with a large uncertainty since a long extrap- 
olation is required. This fraction is also consistent 
with the fraction inferred (Miller & Hamilton 2002) 
from the observed luminous point sources in several 
starburst galaxies. 

Assuming that massive black holes exist in at least 
some GCs, the rate at which white dwarfs pass suf- 
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ficicntly close to a massive black hole to produce ob- 
servable supernovae out to cosmological distances is 
given by a combination of factors, 



NmBHSN = Ng a i X Nqc X Jbh X X\VD 



(15) 



where N ga i is the number of galaxies out tozwl, 
Nqc is the average number of globular clusters per 
galaxy, and Jbh is the fraction of those clusters which 
contain a black hole with tubh > 1000 M . 

For globular-cluster black holes, the rate Xwd at 
which white dwarfs pass sufficiently close to explode 
can be written, 



Xwd = n core (av) 



(16) 



where n core is the density of carbon-oxygen white 
dwarfs around the black hole. The quantity (<rv) 
is the average white-dwarf gravitational cross section 
times the virial velocity distribution in the dense envi- 
ronment near the black hole. Equation (16) assumes 
an isotropic distribution of stars and a relaxation time 
for the system < the rate of black-hole encounters. 

The gravitational cross section a can be written in 
terms of a maximum asymptotic impact parameter b, 



a = irb 2 



(17) 



such that a white dwarf passes close enough to the 
black hole for its density to exceed the pycnonuclear 
ignition threshold. 

First, let us estimate the cross section a. A white 
dwarf infalling from a large distance will have a spe- 
cific angular momentum: 



L = EVryob , 



(18) 



where is the tangential velocity component of the 
white dwarf at a large distance from the black hole. 
The specific mass-energy is given by£~ 1.0+d^/2 . 
This gives 

L = Eb^2(E -1) . (19) 
Now, at closest approach z m in, 



E 



= ^{l-2M/z mm )(l + L 2 /z 2 min ) , (20) 



which can be solved for the impact parameter b to 
give: 



2Mz n 



v 2 00 (l-2M/z min ) ■ 



(21) 



Now inserting, U ~ 8z min = L/z min , leads to 
U 2 w b 2 v 2 lz 2 = — 

oo / ^mtn /i OA/fly \ 



(22) 



For the cases of interest in Tables 2 and 3, Mj ' z m in is 
small at the distance at which the ignition density is 
reached and can be ignored. Solving this for b 2 and 
substituting U 2 w 2M/z min , gives the desired cross 
section as a function of v^, for the white dwarf, M 
for the black hole, and the U 2 at which the star can 
explode, 

,2 4ttM 2 

Equivalently, one can write the cross section in terms 
of the distance of closest approach, 



_ 2irMz min _ 2 / W_ 

u oo \ u oc 



(24) 



Let us consider a typical globular cluster (Binney 
& Tremaine 1987) with a core radius of 1.5 pc, a core 
density of n core w 8 x 10 3 M© pc 3 , and a core virial 
velocity of Voo w 7 km s _1 . Even at high redshift, the 
core is likely to be dominated by white dwarfs with a 
mass distribution peaked around 0.6 M© but with a 
broad distribution extending up to 1.2 M© (Silvestri 
et al. 2001). We have determined the rate at which 
supernovae of this type are ignited per globular clus- 
ter by integrating Eq. (16) over a thermal distribu- 
tion of white dwarf velocities normalized to a virial 
velocity of ^ w 7 km s _1 . The velocity dependent 
gravitational cross section was taken from Eq. (24) 
with an average for z m i n taken from the results of 
Table 2 weighted by the mass distribution of Silvestri 
et al. (2001). We determine in this way a rate of 
Xwd = 2 x 10~ 8 SNe per year per massive black hole 
in a globular cluster. Hence, the use of Eq. (16 is jus- 
tified since our estimated black-hole encounter rate is 
less than the median core relaxation time derived for 
GCs in the Milky Way (Binney & Tremaine 1987). 

The number of galaxies out to a redshift z w 1 
can be obtained by integrating the galaxy luminosity 
function and is approximately N ga i w 10 9 . The num- 
ber of globular clusters per galaxy is a strong function 
of the parent galaxy luminosity (Vandenbergh 1997). 
For example, in M87 there are 1.4 x 10 4 globular clus- 
ters while the Milky Way only contains about 150. 
Since globular clusters are destroyed over time, it is 
quite possible the the low number of GCs in the Milky 
Way is atypical. If we take the number of GCs around 
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M87 as more realistic, then there could be more than 
10 13 globular clusters out to a redshift of z ~ 1. 

As of yet there are no confirmed globular clusters 
which contain a massive black hole in the Milky Way. 
Nevertheless, as argued above it is at least possible 
that a fraction of GC's contain an as yet undiscov- 
ered massive black hole. Therefore, we allow for the 
possibility that Jbh > 0. The whole sky rate for this 
class of supernova then becomes 

Nmbhsn^I^^Jbh yr" 1 . (25) 

Although this estimate (like that for LMBHSNs) is 
smaller than the observed (~ 10 7 yr -1 ) Type la su- 
pernova rate out to large redshift, it could be much 
greater if for example, the average number of globu- 
lar clusters per galaxy were larger in the past. Hence, 
even for Jbh $0.1, we conclude that this class of su- 
pernova event may be detectable. 

5.2. Supermassive Black Holes in Galactic 
Centers 

We next estimate the rate for the explosion of 
white dwarfs on passing orbits near supermassive 
black holes in a galactic centers. Such systems may be 
natural candidates for the events described here. For 
one, it appears that almost any moderate size galaxy 
has a supermassive black hole at its core (Kormendy 
& Richstone 1995; Magorrian et al. 1998) with a mass 
ranging from 10 6 to 10 9 M Q (Gebhardt et al. 2003). 
Furthermore, it is well know that quasars and AGN's 
are absorbing of order 1 to 10 M© of material per year, 
and that there is a very high stellar concentration 
around the central supermassive black hole. For ex- 
ample the modest 3 x 10 6 M© SgrA* black hole in the 
Milky Way has an observed mass concentration after 
subtracting the black hole (Ghez et al. 1998; Genzel 
et al. 2003) of 1.2 x 10 6 x (i?/10") M© pc" 3 where R 
is the angular distance from SgrA* in arcseconds (1 
pc :=y 24" in the Galactic center). In the inner cusp, 
the density can exceed 10 8 M© pc -3 . Although many 
of these stars appear to be members of a young stellar 
population (Figer 2000; Gezari et al. 2002), it seems 
likely that these young stars are the result of mergers 
from older stars (Schodel et al. 2003). Hence, a com- 
parable population of white dwarfs could be present. 
We conservatively estimate n w d ~ 10 6 pc~ 3 around 
the black hole with an average velocity dispersion of 
Woo ~ 100 km s _1 . 



In the case of a galactic center, equation (16) can- 
not be used to estimate the steady state rate at which 
stars fall into the BH. This is because stars on eccen- 
tric orbits that bring them close enough to the BH to 
be destroyed ("loss-cone" orbits) are depleted rapidly 
on a dynamical time scale. Subsequently, the rate is 
set by the much slower diffusion time-scale of stars 
from BH-avoiding orbits into radial orbits. Here we 
will apply the loss cone model as discussed in Syer & 
Ulmer (1999) and Magorrian & Tremaine (1999). 

In this model the two-body relaxation timescale for 
a spherical cluster of stars of mass m* with a density 
p{r) and an isotropic density dispersion a(r) is given 

by 



(26) 



In kG 2 pm* ' 

where m* is the typical stellar mass and A includes 
dimcnsionlcss factors of order unity as well as the 
Coulomb logarithm. Following Binney & Tremaine 
(1987) and Syer & Ulmer (1999) we then can write a 
numerical value of: 



1.8 x 10 10 y 
ln(0.4 x N) V 100 km s" 1 
'M Q \ ^10 4 M© pc~ 3 
rn 



(27) 



where N is the number of stars within the character- 
istic radius Tb and a is the one-dimensional velocity 
dispersion (Seyer & Ulmer 1999) given as 



GM 
(l+p)r ' 



(28) 



where p ~ 2 is the logarithmic slope of the density. 
The white dwarf capture rate will then be given by 
the rate at which white dwarfs will scatter into the 
loss cone. 



A 



WD 



I 



^im wd {r)r 2 dr 



(29) 



In (2/6 lc )t r ' 
where n w d{r) is the number density of white dwarfs 



n w d = fwdp/ (m w d) 



(30) 



where f wc i ~ 0.5 is the mass fraction in white dwarfs 
and (m w d) ~ 0.6M© is the mean white-dwarf mass. 
The quantity 9i c = qGM/r 2 a 2 is the angular size of 
the "loss cone" where q is the radius at which stars 
are removed from the system. The density profile p(r) 
divides into inner and outer regions. The black hole 
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dominates the gravity in the inner region so that a 
density profile of the form p w r~ 2 and a dispersion 
velocity, a 2 w GM/3r is expected. In the outer re- 
gion on the other hand one expects flatter density and 
velocity-dispersion contours. Based upon the analysis 
of Syer & Ulmer (1999) (cf. their figure 4) we esti- 
mate Xwd ~3x 10~ 5 for a typical white dwarf mass 
of 0.6 M©. However, based upon the range and un- 
certainties in the parameters p, a, r^, and Mbh we 
expect a range in Xwd from up to an order of magni- 
tude above to a couple of orders of magnitude below 
this number. For purposes of discussion we will adopt 
~ 3 x as a reasonable estimate. 

The total rate of such events out to z w 1 will be 
given by a similar combination of the factors to those 
given above in Eq. (15), 

RSMBHSN = NgaiNsM^WD , (31) 

where again, N ga i ~ 10 9 is the number of galaxies 
out to z w 1 and Nsm ~ 1 is the number of super- 
massive black holes per galaxy. The quantity Xwd is 
again the rate at which white dwarfs are making suf- 
ficiently close passes to a black hole. The total rate 
is then of order 10 3 to a few xlO 5 yr _1 . This rate is 
somewhat less than normal SNla's which are observed 
at a rate <~ 10 7 yr _1 . Furthermore, a large fraction 
of such events might be hidden in the dense regions 
surrounding galactic cores. Nevertheless, this crude 
estimate is at least suggestive that such events could 
be frequent enough to be observed. It is conceivable 
that they might even dominate over normal Type la 
supernovae at sufficiently large redshift where newly 
formed galactic cores have not yet depleted their ra- 
dial orbits. Frequent white dwarf encounters with a 
supermassive black hole might even provide a power 
source for AGNs and quasars. Clearly, if such explo- 
sions occur, a search for observed signatures seems 
warranted. 

6. Possible Observed Signatures 

Having proposed this class of new supernovae it is 
worth commenting on several distinguishing charac- 
teristics which might help observers identify whether 
one of these events has actually occurred. Prelimi- 
nary numerical burn simulations (D. Dearborn, priv. 
comm.) of 0.6 M white dwarfs passing a supermas- 
sive black hole indicate a nearly complete burn of ini- 
tial C/O into Ni/Co. Hence, each supernova should 
indeed generate ~ 10 51 erg in kinetic energy of mat- 



ter and ~ 10 49 erg in optical light. The optical signal, 
however, may be unique for several reasons. 

One obvious feature is that these supernovae should 
most likely be seen in association with a galactic 
core (or possibly a globular cluster). Indeed, a num- 
ber of transients that might be such supernovae near 
galactic cores have been observed (N. Suntzeff, priv. 
comm.) out to large redshift, but are often not an- 
alyzed due to the complication of their location in a 
crowded field. We suggest that in the near future, the 
development of the National Virtual Observatory (cf. 
Hanisch 2002) might be utilized to search for a pos- 
sible association of unusual Type I supernovae with 
galactic cores and/or extragalactic globular clusters. 

Among other things to look for is the fact that the 
stars are moving very fast when they explode. Hence, 
they could exhibit a large Doppler shift in their spec- 
tra. Also, the explosion itself should be different, 
since the stars are of lighter mass and undergo a very 
rapid rapid compression/ignition. Normal Type la 
events are expected to undergo an off-center ignition 
after a gradual approach to the ignition density. The 
events described here, however, will probably ignite 
rapidly and at their centers. This will cause a det- 
onation and a possibly more efficient thermonuclear 
burn of the white dwarf (depending upon tidal effects 
from the black hole). Thus, one might expect a dif- 
ferent burn front, a different light curve, and different 
nucleosynthesis yields (e.g. more iron and reduced sil- 
icon) compared to normal SNIa events. 

6.1. Interaction with the Black Hole 

Another distinguishing characteristic of these ex- 
plosions is that they occur in the strong field of a com- 
panion black hole. Hence, there might be an affect on 
the light curve from matter accreting into the black 
hole and/or impinging on the accretion disk. For ex- 
ample, the accretion of SN debris onto the black hole 
might produce excess X-ray (or even 7-ray) emission. 

We suggest two possible observational consequences 
of the SMBHSNe, depending upon the optical depth 
of the material in the accretion disk down to the loca- 
tion where the supernova occurs. If the optical depth 
to the point of explosion is greater than one, then 
since the scale of the systems is so large, energy will 
be emitted from the system over a timescale of a year 
or more. 

In the case that the optical depth is less than one 
then the usual SN light signal will be seen. However, 
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if there is an accretion disk present, then the SN light 
signal will be augmented by radiation emitted from 
the collision of the SN remnant with the accretion 
disk. 

Let us consider order-of-magnitude estimates for 
the radiation emitted as debris from the explosion of 
a 0.6 M Q white dwarf impinges on the accretion disk 
of a 10 9 solar mass black hole. For this we will assume 
that the accretion disk is flat. That is, we ignore the 
thickness of the accretion disk near the black hole. 

From Table 3, the distance from the black hole 
to the point at which the explosion initiates is typ- 
ically about d = 1.5 x 10 15 cm (roughly ten black 
hole Schwarzschild radii). We take this as a measure 
of the distance from the supernova to the accretion 
disk. The expansion velocity of the supernova shell 
is given by v w (2 x 10 51 /M W r>) 1/2 = 1-3 x 10 9 cm 
s _1 . Hence, the time scale for material to reach the 
accretion disk is of the order 1.5 x 10 15 /v w 10 6 sec. 
The time for peak luminosity of a Typc-I supernova 
is also about 10 6 sec. Thus, the peak supernova sig- 
nal and the start of the thermal radiation from the 
collision of the supernova remnant with the accretion 
disk will typically occur at about the same time. 

The peak temperature from the collision shock 
will be about 400 keV. The electrons and baryons 
will come to equilibrium temperatures quickly. Then 
bremsstrahlung radiation will cool the hot gas on a 
timescale of about 10 4 sec. Most of the remnant will 
not hit the accretion disk perpendicularly. Averaging 
the various physical quantities over the angle of im- 
pact we estimate a mean photon energy of about 30 
keV. The time for the photon spectrum to fall below 
10 keV will be about 10 7 sec. The total emitted X-ray 
energy will be ~ 10 50 ergs. Thus, an X-ray burst in 
coincidence with a supernova explosion deep inside a 
galaxy is what we expect to possibly be observed. 

6.2. Remnants 

As another observed signature, it is possible that 
remnants from this type of supernova may be de- 
tectable in the Galaxy. For example, an event rate of 
a few xlO -5 y~ 4 per galaxy for white dwarfs to ap- 
proach and explode near the giant black hole in galaxy 
centers, suggests that it is worthwhile to search for a 
supernova remnant with an age w 3 x 10 4 y near the 
giant black hole SgrA* in the Galactic center. For a 
remnant moving with a typical dispersion velocity of 
100 km s _1 for 10 4 y, one would expect to find such 



a remnant <~ 1 pc from SgrA*. 

In this regard, it is certainly of interest that in- 
deed such a remnant (SgrA East) exists and has 
been recently analyzed (Maeda et al. 2002) with 
the ACIS detector aboard the Chandra X-ray Ob- 
servatory. This remnant is centered at a distance of 
w 3 pc from SgrA*. In fact, SgrA* is within the 
remnant outer rim. The passage of the supernova 
shell sometime in the past may have swept gas away 
from the black-hole vicinity and hence, may be re- 
sponsible for the present quiescent state of SgrA*. 
Also of note is the fact that its peculiar metal abun- 
dances, along with the unusual combined radio and 
X-ray morphologies classifies this remnant as a metal- 
rich " mixed- morphology" (MM) SNR. Although, the 
properties of SgrA East may be explained (Maeda 
ct al. 2001) by a low-mass Type II supernova, it 
is tempting here to speculate that this class of MM 
SNRs may be the smoking gun for a black-hole in- 
duced white-dwarf explosion (cf. Khokhlov & Melia 
1996). 

Obviously, a clear picture of how to classify these 
events will require fully relativistic numerical compu- 
tations in three spatial dimensions of the compres- 
sion and thermonuclear explosion of a white dwarf in 
the background field of a black hole. Efforts along 
this line are currently under way (D. Dearborn priv. 
comm.) utilizing three dimensional hydrodynamics 
and radiation transport. As noted above, preliminary 
results indicate a central thermonuclear ignition (very 
close to our adopted ignition density) which proceeds 
to detonate the star. In future work we will simulate 
the radiation transport, light curve, and supernova 
remnant. 

7. Conclusions 

We have explored a new class of Type-I supernovae 
whereby the onset of a thermonuclear explosion is in- 
duced by relativistic enhancements of the white-dwarf 
self gravity as it accelerates in the background grav- 
itational field of a black hole. We find that a po- 
tentially observable rate of such events could be oc- 
curring out to cosmological distances, particularly in 
the dense regions of galactic cores and possibly glob- 
ular clusters as well. These explosions and their rem- 
nants could thus be characterized by their association 
with supermassive black holes galactic cores or mas- 
sive black holes in globular clusters. They also might 
be characterized by an associated x-ray burst and/or 
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a significantly modified light curve from matter col- 
liding with the black-hole accretion disk. 
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Fig. 1. — Contours of constant density in the r vs. z 
plane for a 0.6 Mq white dwarf approaching a 1000 
M Q black hole. The star is near closest approach, 
centered at z 1.4 x 10 4 km from the black hole 
which is located below the figure. This star is tidally 
distorted to an aspect ratio of (z, r) w (900, 500) km. 
Nevertheless, the radius has contracted by a factor of 
10 from that of an isolated star. The central density 
has correspondingly increased by more than a factor 
of 1000 from (3.2 x 10 6 ) to (3.3 x 10 9 g cm" 3 ) which 
is above the pycnonuclear ignition threshold. 
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Fig. 2. — Contours of constant density in the r vs. z 
plane for a 0.6 M© white dwarf near closest approach 
(z k, 2 x 10 4 km) to a 1500 M black hole. 



Fig. 3. — Contours of constant density in the r vs. z 
plane for a 0.6 M white dwarf near closest approach 
(z«4x 10 4 km) to a 3000 M Q black hole. 
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